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Abstract 

The energy density e g of asymptoticaly flat gravitational fields can be calculated 
from a simple expression involving the trace of the torsion tensor. The integral of 
this energy density over the whole space yields the ADM energy. Such energy 
expression can be justified within the framework of the teleparallel equivalent of 
general relativity which is an alternative geometrical formulation of Einstein's gen- 
eral relativity. In this letter we apply e g to the evaluation of the energy per unit 
length of a class of conical defects of topological nature, which include disclinations 
and dislocations (in the terminology of crystallography). Disclinations correspond 
to cosmic strings, and for a spacetime endowed with only such a defect we obtain 
precisely the well known expression of energy per unit length. However for a pure 
spacetime dislocation the total gravitational energy is zero. 
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I. Introduction 

It is believed that phase transitions in the early universe can give rise to topological defects, 
which can lead to very important cosmological consequences]]]]. In order to understand 
the formation of galaxies and cluster of galaxies it has been suggested that these structures 
have evolved from the gravitational instability of small density fluctuations. One of the 
major unresolved cosmological problems is the origin of these initial density fluctuations. 
One possibility is that the latter are due to cosmic strings 0, [|. For this reason cosmic 
strings have been widely studied in the literature. 

In the terminology of crystallography^, [j| cosmic strings correspond to disclinations, 
which is one possible defect in a crystal. To our knowledge, other commom crystal defects 
like dislocations have not yet been considered in the various cosmological models. From 
a geometrical point of view, disclinations and dislocations are conical singularities in 
a flat, four dimensional Lorentzian spacetime, i.e., they can be described by a metric 
which is flat away from the r = axis, but with a coordinate singularity that cannot be 
removed. Tod0 has recently considered these defects in a unified fashion. He generalized 
an argument due to Vickers0, according to which a conical singularity like a cosmic 
string can be interpreted in terms of a ^-function of curvature supported in the origin of 
the (x, y) plane, say. Tod extended this idea to conical sigularities of the dislocation type 
and argued likewise that instead of interpreting the latter as defects in an otherwise flat 
Minkowski spacetime, one can alternatively consider a flat Minkowski spacetime endowed 
with a delta function of torsion as a source in the r = axis. 

In this paper we will calculate the gravitational energy of the field configuration con- 
sidered by Tod, which includes altogether disclinations and dislocations. We will obtain 
the energy per unit length along the defect axis r = 0. For a metric field which describes a 
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disclination only the energy per unit length turns out to be exactly the same of the cosmic 
string. Together with other previous calculations of energy of black hole configurations, 
this result supports the validity of the present energy expression. 

For the dislocations the result is remarkable. The energy per unit length of these 
defects is zero, when the integration extends over the whole three dimensional space. 
This result may be of significance for cosmology. It may be conjectured that a physical, 
cosmological discolation needs a very low energy to be formed, and hence should play an 
important role in phase transitions in the early universe. 

The difficulty in obtaining an expression for localized energy density in the framework 
of the Hilbert-Einstein action integral has led to the widespread belief that the gravi- 
tational energy cannot be localized. We do not share this idea. The energy density of 
the gravitational field can be naturally obtained from the Hamiltonian formulation of the 
teleparallel equivalent of general relativity (TEGR) || . The TEGR is an alternative geo- 
metrical formulation of Einstein's general relativity. The gravitational field in the TEGR 
is described by the tetrad field, and its dynamics is dictated by Einstein's equations. 
Therefore this is not an alternative theory of general relativity. The gravitational energy 
density for asymptotically flat geometries has been presented and justified in ref.0. 

In section II we present the mathematical preliminaries of the TEGR, its Hamilto- 
nian formulation and the expression of the energy for an arbitrary asymptotically flat 
spacetime. We also calculate the energy inside a surface of constant radius r Q for both 
the Schwarzschild and the Kerr solutions. In section III we present the calculation of the 
energy of conical defects. 



II. The TEGR in Hamiltonian form 

Notation: spacetime indices and local Lorentz indices a,b,... run from to 3. 
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In the 3+1 decomposition latin indices from the middle of the alphabet indicate space 
indices according to \i = 0,i, a = (0), (i). The tetrad field e a ^ and the spin con- 
nection u^ ab yield the usual definitions of the torsion and curvature tensors: R a bfJiU = 
d fJl u> u a b + oj^ a c oj u c b — T a MI/ = d^e a v + be b v — .... The flat spacetime metric is 
fixed by r\ {m) = -1. 

In the TEGR the tetrad field e a M and the spin connection u^ ab are completely inde- 
pendent field variables. The latter is enforced to satisfy the condition of zero curvature. 
The Lagrangian density in empty spacetime is given by[Q] 

L{e,to,X) = -ke(- A T abc T abc + ]f abc T bac - T a T a ) + e\ ab ^R ab ^{u) . (1) 

where k = j^q, G is the gravitational constant; e = det(e a M ), \ ab ^ u are Lagrange 
multipliers and T a is the trace of the torsion tensor defined by T a = T b ba . 

The equivalence of the TEGR with Einstein's general relativity is based on the identity 

eR{e,u) = eR(e) + e{-T abc T abc + l -T abc T acb - T a T a ) - 2d,(eT») , (2) 

which is obtained by just substituting the arbitrary spin connection uj^ ab = °0J lJiab (e) + 
K^ ab in the scalar curvature tensor R(e,uj) in the left hand side; °uj^ ab (e) is the Levi- 
Civita connection and K^ ab — |e x e b u {T\ fJiU + T u \^ — T^x) is the contorsion tensor. The 
vanishing of R a b/iu (u), which is one of the field equations derived from (1), implies the 
equivalence of the scalar curvature R(e), constructed out of e% only, and the quadratic 
combination of the torsion tensor. It also ensures that the field equation arising from the 
variation of L with respect to e a M is strictly equivalent to Einstein's equations in tetrad 
form. Let -M^ = denote the field equation satisfied by e a „. It can be shown by explicit 
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calculations that 



SL 1 r 1 

:{Ran - -e a ^R{e)\ . 



5e a» 2 L M 2 
(we refer the reader to refs.|| |9| for additional details). 

It is important to note that for asymptoticaly flat spacetimes the total divergence 
in (2) does not contribute to the action integral. Therefore the latter does not require 
additional surface terms, as it is invariant under coordinate transformations that preserve 



the asymptotic structure of the field quantities ||T0| . 

The Hamiltonian formulation of the TEGR can be successfully implemented if we 
fix the gauge u 0a b = from the outset, since in this case the constraints (to be shown 
below) constitute a first class set||. The condition a> 0a 6 = is achieved by breaking 
the local Lorentz symmetry of (1). We still make use of the residual time independent 
gauge symmetry to fix the usual time gauge condition e^) = e( )j = 0. Because of 
^oab = 0, H does not depend on p kab y the momentum canonically conjugated to u>kab- 
Therefore arbitrary variations of L = pq — H with respect to P kab yields iokab — 0. Thus 
in view of c^oafe = 0, u^ab drops out from our considerations. The above gauge fixing can 
be understood as the fixation of a global reference frame. 

Under the above gauge fixing the canonical action integral obtained from (1) becomes|3[] 

A TL = J d A x{H^ k e m - H} , (3) 

H = NC + N*Ci + S mn n m " + J—d k (NeT k ) + ^(IF'*^) . (4) 

N and are the lapse and shift functions, H mn = eu\ m Uy) n and S mn = — S nm are 
Lagrange multipliers. The constraints are defined by 
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C = d^keV) - keE k »T kij - — (H'-MI,; - -II 2 ) , (5) 

C k = -e u)k d t U^ - U^% )ik , (6) 

with e = det(e(j) k ) and T % = # lfc e (j) 'T (j)Zfc , T (j)Zfc = die^ j)k - <9 fc e (i)i . We remark that (3) 
and (4) are invariant under global S0(3) and general coordinate transformations. 

We assume in this section the asymptotic behaviour e^) k ~ r)j k + \hj k {\) for r — > oo. 
In view of the relation 

( (PxdjieT*) = f dSkidihik-dkha) = E ADM (7) 

87rG J lbnG Js 

where the surface integral is evaluated for r —>■ oo, we note that the integral form of the 
Hamiltonian constraint C = may be rewritten as 

J d*x{keY, Hj T m + -^(ir'll,, - ±n 2 )} = E ADM . (8) 

The integration is over the whole three dimensional space. Given that <9j(eT J ) is a scalar 
density, from (7) and (8) we define the gravitational energy density enclosed by a volume 
V of the space as[|9] 

e b = -Trf d 3 xdj(eT3) . (9) 
87rCx Jv 

It must be noted that this expression is also invariant under global SO (3) transformations. 

We will briefly recall two applications of E g . Let us initially consider a spherically 
symmetric geometry and fix the triads ea.v as 
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/ e A sinOcoscj) r cosOcoscp —r sinOsincf) \ 

£(k)i — e A sinOsincf) r cosOsincp r sinOcoscp 

\ e A cosO — r sinO 
(k) is the line index and i is the column index. The function A(r) is determined by 



(10) 



e- 2A = 1 



2mG 



The one form e( fc ) = r dr + gd6 + ^dcfr yields 



e (k) e {k) 



e 2A dr 2 + r 2 de 2 + r 2 sm 2 6i 



Therefore the triads given by (10) represent the spatial section of the Schwarzschild solu- 
tion. We can easily calculate e g = ■^^d i {eT l ) associated to (10). We obtain 



]_d_ 



r(l - e- A )] , 



The energy inside a spherical surface of arbitrary radius r Q can be calculated from (11). 
It is given by|§ 



E n 



rjl 



1 



2mG, i 



(12) 



This is exactly the expression found by Brown and York|TT[ in their analysis of quasi-local 
gravitational energy. They define a general expression for quasi-local energy as minus the 
proper time rate of change of the Hilbert-Einstein action (with surface terms included), 
in analogy with the classical Hamilton- Jacobi equation which expresses the energy of a 
classical solution as minus the time rate of change of the action. The application of their 
procedure to the Schwarzschild solution yields (12). Note that when r Q — >• oo we find 
E = m. 



The definition (9) for the gravitational energy can also be successfully applied to the 
Kerr black hole[T2[]. In terms of Boyer and Lindquist coordinates |TB[ (t, r, 9, 4>) the spatial 
section of the Kerr metric is given by 



ds 2 



with the following definitions: 



n 2 Y? 

—dr + p dO H -sin 

A p 2 



(13) 



A 



r 2 — 2mr + a 2 , 



2 2,2 2/1 

r = r + a cos v , 



( r 2 + a 2)2 _ Aa 2 sm 2^ 



a is the specific angular momentum defined by a = — . The triads appropriate to the 
three-metric above are given by 



/ -^sinO coscj) pcosO cos<p — ^sin8 sin<j)\ 



-^=sin8 sintfi pcosO sincj) ^sinOcosc, 



/ 



(14) 

we easily find 



—psind 
Indeed, defining again the one-form = e^ k ' r dr + e^ k ' gd0 + 

that e^e(fc) = (is 2 given by (13). 

There is another set of triads that yields the Kerr solution, namely, the set which is 
diagonal and whose entries are given by the squere roots of gu. This set is not appropriate 
for our purposes, and the reason can be understood even in the simple case of flat space- 
time. In the limit when both a and m go to zero (14) describes flat space: the curvature 
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and the torsion tensor vanish in this case. However, for the diagonal set of triads (again 
requiring the vanishing of a and m), 



M = dr , e {e) =rd9 , e w = r sin9 . 



some components of the torsion tensor do not vanish, Tm)%2 = 1, Tt 3 )i 3 = sin9, and E g 
calculated out of the diagonal set above diverges when integrated over the whole space. 
Moreover, it is the flat space form of (14), i.e., when m = a = 0, that can be brought 
to a diagonal form in cartesian coordinates, and not the diagonal form above in spherical 
coordinates. Thus the asymptotic behaviour ~ rjjk + \h{-) when r — > oo can only 
be achieved by means of (14). 

In ref. [II] we have obtained the expression of the energy contained within a surface of 
constant radius r = r D : 



E„ 



d6 sin6< p 



2r(r 



2 ' a 2 ) 



a 2 sin 2 6(r — m) 



(15) 



4 Jo ' V ' p pS 

In the limit of slow rotation, namely, when ^ << 1 all integrals in the expression 
above can be calculated, and E g finally reads 



E n 



Toil 



\ 



2m a 2 " 

1 + ^7 



\ ° 2 




2m 






2H 


+ 




/ 6r 









2m a 2 " 
1 + - 



(16) 



This is exactly the expression found by Martinez who approached the same problem 
by means of Brown and York's procedure. However the present approach is more general 



than that of ref.[H)|. The energy given by (15) can be calculated by means of numerical 
integration for any value of a. On the other hand Brown and York's procedure requires 
the embedding of an arbitrary two dimensional surface of the Kerr type in the reference 
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space E 3 , a construction which is not possible in general (T^j (the evaluation of the energy 
in ref.[I5| is only possible in the limit — << 1). 



III. Conical Defects 

The calculations of the previous section support the correctness of expression (9) for 
the energy of the gravitational field. In ref.§J expression (9) was justified in the framework 
of asymptoticaly flat gravitational fields. The action integral of the TEGR for compact 
spacetimes differs from the one for asymptoticaly flat geometries by a surface (bound- 
ary) term and consequently the two Hamiltonian densities also differ by a surface term. 
However the Hamiltonian constraint is the same for both kinds of geometries (except, 
of course, for the possible presence of additional terms; had we added a cosmological 
constant to the Lagrangian density, such a term would also appear in the Hamiltonian 
constraint). 

We have seen from (8) that the integral form of the Hamiltonian constraint equation 
may be written as C = H — Eadm = for asymptoticaly flat spacetimes. We will assume 
here that this is a general feature of gravitational theories, namely, we will assume that 
for an arbitrary geometry the Hamiltonian constraint equation may be written in integral 
form as C = H — E = 0. Therefore we will tentatively evaluate the energy of simple well 
known geometries which are not asymptoticaly flat by means of (9). In the following we 
will evaluate the latter for the class of geometries considered by Tod||. Such a class of 
geometries is particularly suitable for our purposes since the energy per unit length of a 
cosmic string is a priori known. 

We will consider conical singularities along the z axis, in an otherwise flat Minkowski 
spacetime, described by the metric 
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ds 2 



-(dt + ad(f>y + dr 2 + (3 2 r 2 d(b 2 + (dz + 7^) 



(17) 



where a, (3 and 7 are real constants and <fi runs from to 2n. The metric is everywhere 
flat except at the axis r = 0. For a = 7 = and (3 < 1 the metric describes a cosmic 
string type singularity. Thus (3 parametrizes a disclination. As shown in ref.|J, a and 7 
parametrize dislocations, in the terminology used in crystallography. The metric as given 



above has also been considered by Gal'tsov and Letelier|T6[. The spatial section of (17) 
reads 



9ij 



(I 0\ 

5 2 7 



VO 7 1/ 

where <5 2 = /? 2 r 2 + 7 2 — a 2 . The corresponding triads are given by 



/ coscj) —5sincb —^siri(j)\ 



e( k )j 



sine 



5 cos ( 



:COSL 



(19) 



V V 1 -^ 

Recall that (k) and j are the line and column index, respectively. Note that if we 
make a = 7 = 0,/?=l(19) can be brought to a diagonal form in cartesian coordinates 
by means of a coordinate transformation. Initially we evaluate the components of the 
torsion tensor: 



12 



or 



— [coscb] = (1 — S^sinc, 



Tans = - 



7 • 
-~sin(i 





— d sine. 
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r (1)23 = ^{- r int) - Q- z {- 5sin v = -r os(t) 



- (2)12 



— (1 — 5') cose 



P (2)13 = -—SC0S(, 



5 2 



7 

^(2)23 = --^sinc, 



T(3)12 = 



O 12 



f - <5 2 

J r3H3 = — ; ^ 



A 3 



^(3)23 = , 

where the prime denotes differentiation with respect to r. We wish to evaluate (9) and for 
this purpose we need the expression of T\ After a long but simple calculation we obtain 

T 2 = T 3 = . 

Together with e = V^ 2 — 7 2 , the energy density can now be easily obtained and integrated. 
The energy E g contained within a cylindrical region with length L and radius r is given 
by 
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We will consider next the three individual situations in which the metrics are parametrized 
by only one of the parameters. 



I. a = 7 = 

The metric parametrized by j3 only describes a disclination. Expression (20) reduces to 

E g = ±{l-(3) (21). 

This is precisely the energy per unit lenght for a cosmic string 0. We note that E g above 
does not depend on the radius of integration r Q . Therefore we may conclude that the 
whole energy is concentrated along the defect axis r = 0. 

II. a = 0, P= 1 

In this case we have a simple dislocation parametrized by 7. We are mostly interested in 
the value of E g for very large values of r Q . From (20) we obtain 

E ^-^h- (22) 

Therefore in the limit when both L and r go to infinity, namely, when the integration is 
performed over the whole three dimensional space, E g vanishes. Thus the total energy of 
in the disclocation is zero. However, in the limit r Q — > we find E g — > —h. 

III. /3 = 1,7 = 
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For this metric expression (20) reduces to 



2 ■ 



For large values of r Q we obtain 

^»-t£ < 23 >- 

Therefore the total energy corresponding to this disclination is also zero. In the limit 
r„^awe observe that E g — > — oo. 

We observe that whereas for a disclination the whole energy per unit length is concen- 
trated along the defect axis, for both types of dislocations the energy is distributed over 
the whole three dimensional space. 



IV. Comments 

By means of expression (9) for the energy of the gravitational field we have obtained 
the correct value of the energy per unit length of a cosmic string. This fact is a good 
indication that the Hamiltonian constraint may be generically written as C = H — E = 0, 
which is the assumption we made in this paper. It also supports the conjecture that 
expression (9) might have a universal character, since it yields the expected values of 
energy for totally distinct spacetimes, namely, Kerr-type and conical spacetimes. 

From the results of section III we conclude that dislocations are more likely to appear 
as a result of cosmological phase transitions in the early universe than disclinations (cosmic 
strings). The energy of an actual, physical dislocation might be nonvanishing, but anyhow 
we would expect it to be very small. The situation here is very much similar to what 
happens in a real crystal. It is well known that disclinations require too much energy to 
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be formed in crystals, whereas dislocations are more favourable defects since they require 
much less energy (see, for instance, sections 6.3.2 and 6.5 of ref.[[5J for a discussion as 
to why the energy cost for a disclination in crystal is prohibitively high). In order to 
understand the vanishing of the total gravitational energy for dislocations, let us consider 
the metric for which a = 0,(3 = 1 (case II). As discussed by Tod||, this metric can be 
transformed into a flat metric if we define a new coordinate Z = z + Z Q ^ = z + 70. 
The associated Burgers vector is determined by Z = 2ttj. We know, however, that the 
energy of an actual dislocation in a crystal depends not only on the Burgers vector but 
also on an extra quantity, the rigidity modulus fi. Such quantity is not given in (17). This 
fact might explain why the total energy of the dislocations above is zero. Disclinations 
and dislocations are concepts used in the deformation models of crystals and metals. We 
conclude from our analysis that dislocations of type considered here might be as well 
useful concepts in cosmological models. 
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